Astronomy & Astrophysics manuscript no. 4571 


February 5, 20l08 


(DOI: will be inserted by hand later) 





o 

> 



Laboratoire d'Astrophysique de Marseille, UMR 6110 CNRS-Universite de Provence, BPS, 13376 Marseille Cedex 12, 
France 



The VIMOS-VLT Deep Survey 

\o'- "^^^ evolution of galaxy clustering per spectral type to z 1.5 * 

o : 

^ ■ B. Meneux' , O. Le Fevre' , L. Guzzo^, A. PoIIqI ^, A. Cappi^ O. Ilbert'^, A. Iovino^ C. Marinoni^ '', H.J. 

McCracken"^ ^ D. Bottini'', B. Garilli'^, V. Le Brun^ D. Maccagni'', J.P. Picat^°, R. Scaramella", M. Scodeggio'^, L. 
^ : Tresse\ G. Vettolani^^^ ^ Zanichelli^^^ q Adami\ S. Amouts\ M. Amaboldi^^ S. Bardelli^ M. Bolzonella"^, S. 
lJ_^ Charlot''^"^, P. Ciliegi^ T. Continii", S. Foucaud^^ P Franzetti'', I. Gavignaudi"-'^, B. Marano"^, A. Mazure\ R. 

Merighi^ S. Paltani^^''^ R. Pello"', L. Pozzetti^ M. Radovich^^ G. Zamorani^ E. Zucca^ M. Bondi^^^ ^ 
^ ■ Bongiorno**, G. Busarello'^, O. Cucciati^'^'^, L. Gregorinii^, R Lamareille^", G. Mathez^°, Y. Mellier^-^, P. 
I Merluzzi ^ ^ , V. Ripepi ' ^ , D . Rizzo ' ° 

(N 
> 

, INAF-Osservatorio Astronomico di Brera - Via Bianchi 46, 1-23807, Merate, Italy 

^ INAF-Osservatorio Astronomico di Bologna - Via Ranzani,l, 1-40127, Bologna, Italy 
Universita di Bologna, Dipartimento di Astronomia - Via Ranzani, 1, 1-40127, Bologna, Italy 
VO . ' INAF-Osservatorio Astronomico di Brera - Via Brera 28, Milan, Italy 

' * Centre de Physique Theorique, UMR 6207 CNRS-Universite de Provence, case 907 F- 13288 Marseille, France 

(-H ^ ' Institut d'Astrophysique de Paris, Universite Pierre et Marie Curie, UMR 7095, 98 bis Bvd Arago, 75014, Paris, France 

* Observatoire de Paris, LERMA, 61 Avenue de I'Observatoire, 75014 Paris, France 
' INAF - lASF Milano, via Bassini 15, 1-20133 Milano, Italy 
O _ Laboratoire d'Astrophysique de I'Observatoire Midi-Pyrenees (UMR 5572) - 14, avenue E. Belin, F31400 Toulouse, France 

" INAF-Osservatorio Astronomico di Roma - Via di Frascati 33, 1-00040, Monte Porzio Catone, Italy 
'2 IRA-INAF - Via Gobetti,101, 1-40129, Bologna, Italy 

INAF-Osservatorio Astronomico di Capodimonte - Via Moiariello 16, 1-80131, Napoli, Italy 
''^ Max Planck Institut fur Astrophysik, 85741, Garching, Germany 
k>( ^ " School of Physics & Astronomy, University of Nottingham, Nottingham NG7 2RD, England 

j_j . European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748, Garching bei Munchen, Germany 

' " Integral Science Data Centre, ch. d'Ecogia 16, CH-1290 Versoix, Switzerland 

Geneva Observatory, ch. des Maillettes 51, CH-1290 Sauverny, Switzerland 
" Universita di Milano-Bicocca, Dipartimento di Fisica - Piazza delle Scienze, 3, 1-20126 Milano, Italy 
email: baptiste.meneux@oamp.fr 



Received ... 2005 ; accepted ... 2005 



Abstract. We measure the evolution of clustering for galaxies with different spectral types from 6495 galaxies with 17.5 < 
Iab ^ 24 and measured spectroscopic redshifts in the first epoch VIMOS-VUT Deep Survey (VVDS). We divide our sample 
into four classes, based on the fit of well-defined galaxy spectral energy distributions on observed multi-color data. We measure 
the projected correlation function Wpir,,) and estimate the best-fit parameters for a power-law real-space correlation function 
^(r) = (r/r^y'. We find the clustering of early-spectral-type galaxies to be markedly stronger than that of late-type galaxies at 
all redshifts up to z 1.2. At z ~ 0.8, early-type galaxies display a correlation length ro = 4.8±0.9 h^'Mpc, while late types 
have ro = 2.5±0.4 h^'Mpc. For the latest class of star-forming blue galaxies, we are able to push our clustering measurement 
to an effective redshift z ~ 1.4, for luminous galaxies (Mb(AB) =^ -21). The clustering of these objects increases up to 
rg = 3.42±0.7 h^'Mpc for z = [1.2, 2.0]. The relative bias between early- and late-type galaxies within our magnitude-limited 
survey remains approximately constant with b = 1 .6±0.3 from z = to z = 1.2. This result is in agrement with the local findings 
and fairly robust against different way of classifying red and blue galaxies. When compared to the expected linear growth of 
mass fluctuations, a natural interpretation of these observations is that: (a) the assembly of massive early type galaxies is already 
mostly complete in the densest dark matter halos at z - 1 ; (b) luminous late-type galaxies are located in higher-density, more 
clustered regions of the Universe at z - 1.5 than their local low luminous counterpart, indicating that star formation activity is 
progressively increasing, going back in time, in the higher-density peaks that today are mostly dominated by old galaxies. 
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1. Introduction 

The measurement of the evolution of the clustering of galax- 
ies is one the key elements to understand the evolution of 
the Universe and the formation of the large-scale structures of 
the Universe. In the current paradigm, the structures we ob- 
serve today in the large-scale distribution of galaxies are the 
result of the gravitational amplification of primordial overden- 
sities in the dark matter density field. Smaller scales are sup- 
posed to collapse first, with virialized structures as galaxies, 
groups and clusters subsequently built up via multiple merg- 
ing, at a rate of growth essentially governed by the value 
of the mean matter density Q^. In this hierarchical forma- 
tion scenario, the evolution of the clustering o f dark matter 
halos is now well understood both analytically dMo & White! 
Il996t ISheth& Torrnenlll999l) and fro r n N-body simulations 
/Jenkins et all Il998t iKauff^mann et all Il999al ISpringel et all 
boQS) . 

The visible component of galaxies is expected to assemble 
via the collapse of baryonic matter within the dark matter halos. 
A naive picture would let to conclude that galaxies follow the 
underlying dark matter density field. The formation and evolu- 
tion of galaxies is, however, driven by physical processes such 
as cooling, star formation and feedback, which are not easily 
linked to the underlying dark matter halos JlCauffman n et al.l 
Il999ah . Rather, a general prediction of most current models 
is that galaxy formation is favored in higher-density regions. 
Objects forming in regions of larger-than-average density col- 
lapse first and, looking back in time, are found to inhabit higher 
and higher density peaks (i.e., they are more biased). Recent 
direct measure ments of the evolution of galaxy bias up to a 
redshift z ~ 1.5 (Marinoni et al. 2005) provide supporting evi- 
dence to this biased picture of galaxy formation. 

Although being well known since more that 25 years, the 
origin of the relationship between galaxy morphological (or 
almost equivalently, spectr al) type and local structure is still 
not fully understood (e.g. IPressleillTgil LSmith et alJEool 
IPostman et al.l EoOSi) . Whether the higher fraction of early- 
type galaxies in high- density regions is established ab ini- 
tio (the so-called nature hypothesis), or rather produced by 
environmental eff'ects during the life of the galaxy (nurture) 
is still debated. A practical effect of the morphology-density 
relation is that early- and late-type galaxies trace differently 
the density field, with the former class being more clustered 
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(i.e. more biased) than the latter. In addition, more luminous 
galaxies (which are preferentially early-types) are more clus- 
tered than less luminous ones, with the two parameters - 
luminosit y and spectral type - playing together in a subtle 
way dlovino et alJll9 93: Mad gwick et alJl20 03: Norbergetal. 
120021 feehavi et al.' 2002; Lovedav T resse & Maddox. .1999i: 
Guzzo et al. 1997; Benoist et al. 1996). Since the amount of bi- 
asing should be related to the mass and formation history of a 
galaxy, one therefore expects the evolution of the clustering of 
galaxies with different morphological types to contain impor- 
tant clues on the assembly of the baryonic mass and the origin 
of galaxy types. 

Large deep surveys consistently show that the clustering 
of the overall galaxy population evolv es slowly with time 
dCoil et alJl2004 iLe Fevre et aliboolt ). Measurements from 
the VIMOS-VLT Deep Survey rVVDS. iLe Fevre et alj2005bl) . 
indicate that the clustering amplitude of L* galaxies increases 
by a factor ~ 2.4 from z - 1.5 to the present epoch. However, 
little is known about the clustering evolution of different galaxy 
types above z ~ 0.3. In the CFRS sample, no significant dif- 
ference was observ ed between the clus tering of red and blue 
galaxies up to z ~ 1 jLe Fevre et alJl996 ). although the sample 
was too small to reach firm conclusions. Red galaxies at z ^ 0.6 
in the CNOC2 survey were found to ha ve a clustering ampli - 
tude ~ 5 times higher than blue galaxies dCarlberg et alJl997l) . 
Results from the CADIS survey using photometric redshifts in- 
dicate that up to z ~ 1 early spectral type galaxies are more 
clustered, with a correlation amplitude 2. 1 times larger than late 
type at z=[0.75-1.07] (Phleps & Meisenheimer 2003). A more 
recent and accurate estimate at (z) = 0.6 ( Phleps et alJl2005l) 
confirms the same trend among red and blue galaxies with 
r'f^''' = 5.39+°:^° h-'Mpc and r^'"' = 3.64+°;25 h-'Mpc. At even 
higher redshifts, the measurement of type is substituted by the 
selection of specific sub-populations, e.g. iDaddi et'al] ( l2003h 
measure a strong clustering of extremely red (J-K) galaxies at 
2 < z < 4, arguing that these may represent the progenitors of 
massive elliptical galaxies at lower redshifts. 

In this paper we make an attempt to quantify accurately 
the contribution of different galaxy types to the observed evo- 
lution of clustering using the VVDS first epoch survey. The 
accuracy of spectroscopic redshifts is used to project the corre- 
lation function ^(rp,7i) in narrow velocity intervals and to mea- 
sure the correlation length ro and the slope y of the real space 
correlation function as a function of redshift up to z ~ 1.2. 

We describe the sample and our galaxy classification in 
Section 2, and the method used to compute the real space cor- 
relation function in Section 3. The clustering measurements 
are presented in Section 4, and then discussed in Section 5. 
We use a cosmology with - 0.3, Qa = 0.7, through- 
out the paper. All distances, given in comoving units, and ab- 
solute magnitudes are computed with h - I (which denotes 
Ho/lOOkms-' Mpc"'). 
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2. VVDS first epoch data 

2.1. The sample 

We are using the VIMOS-VLT D eep Survey first epoch data 
described in iLe Fevre et alJ ( 1200 5 ai. We concentrate on the 
VVDS-02h field, and we keep all galaxies with the best red- 
shifts available, with quality flags 2 to 9 (confidence level 
greater than 80% in the redshift measurement). AGN are ex- 
cluded from the sample. The analyzed sample in this paper in- 
cludes 6495 galaxies with spectroscopic redshifts in an area 
of 1750 arcmin^. The accuracy of the redshift measurement is 
~ 275 km.s''. A detai led description of the sample is given in 
iLe Fevre et all ( l2005ah . 



2.2. Galaxy classification 

We us ed the galaxy classification introduced by IZucca et'aP 
(I2OO5I) to study the evolution of the luminosity function per 
galaxy type from the first epoch VVDS data. The classifica- 
tion has been obtained from the B to I band m ulti-wavelength 
dataset available in the Vy DS-02h field (iLe Fevre et alJ 
|2004aVMcCracken et al.'2003'), and the spectroscopic redshifts 
(Le Fevre et al. 2005a). The U band, which does not cover all 
die field jRadovich et aljboosl) was not used in order to get 
a homogeneous classification on the sky. Rest-frame colors 
have been matched to the empirical set of spectral ene r gy dis - 
tribution (SED) templates describe d in lArnouts et alJ ll 19991) 
made of four observed spectra ("Coleman Wu & Weedman' 
[ l980 , C WW hereafter) and two sta rburst SEDs computed with 
GISSEL llBruzual & Charloil993l) . 

The original spectroscopic sample of 6495 galaxies was 
divided into four spectral classes from early-type (type 1) to 
irregular / star-forming galaxies (type 4) based on the CWW 
templates. The latest type also includes galaxies with SEDs 
well represented by starburst templates. The properties of these 
four classes of galaxies are summarized in Table and their 
redshift distributions are shown in Figure As a result of the 
VVDS I-band selection, we are increasingly loosing red early- 
type galaxies for redshifts z>~ 1 .2 (Figure|3 because galaxies 
are selected from their rest frame UV flux at these redshifts, 
but late type galaxies continue to be identified up to the highest 
redshifts discussed in this paper. 

As described in IZucca et'aP ( l2005l) . the evolution of the 
global luminosity function in the VVDS is strongly driven 
by type 4 galaxies, which evolve by ^ -2 magnitudes from 
z ~ 1 .5 and present a steep faint end luminosity function slope. 
Conversely, the luminosity function of early type galaxies re- 
mains stable over this redshift range. The relative distribution 
of absolute magnitudes per type presented in Figure |3l is the 
result of this evolution. 

In summary, the type classification used in this paper is a 
spectral classification based on a set of spectral templates cov- 
ering from early type galaxies to late and star-forming galaxies. 
Although it is tempting to associate spectral type to morpho- 
logical type, and a close relationship exists between them, the 
spectral classification used in this paper is more closely related 
to a color classification from "red" to "blue" galaxies. 
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redshift z 

Fig. 1. Redshift distribution of the four galaxy types 




0.5 1 1.5 

(B-I)jg (rest-frame) 

Fig. 2. Rest-frame color (B - I)ab for early-type galaxies (type 
1&2) and late-type galaxies (type 3&4) computed from the best 
fit of our multi-wavelength datase t with PEGASE.2 templates 
dFioc & Rocca-Volmerangell997l). 

3. Measurement and associated errors: method 

The method to recover the correlation length ro and the slope 7 
of the real- space correla tion function is well detailed in the pa- 
per of Polio et al. ('2005V We are focusing here on the most crit- 
ical aspects of the method when applied to compute the clus- 
tering properties as a function of galaxy type. 



4 



B. Meneux et al.: VVDS: Clustering evolution per galaxy type 



Table 1. The four VVDS galaxy types 



type 




number of 
objects 


indicated 
rest-frame color 


1 


E/SO 


645 


1.3 < (Bab- Iab) 


2 


early spiral 


1004 


0.95 < (Bab -Iab)< 13 


3 


late spiral 


2104 


0.68 < (Bab - Iab) < 0.95 


4 


irregular 


2742 


(Bab - Iab) < 0.68 




starburst 








Fig. 3. B-band luminosity distribution Mb(AB) for early-type 
galaxies (types 1&2) and late-type galaxies (types 3&4) 

3.1. The real space correlation functions parameters 

We computed the bi-dimensional two-point correlation func- 
tion ^(rp,ji) for each galaxy type, in different redshift bins cho- 
sen to maximize the number of objects hence the signal to noise 
ratio of the clustering mea surement. This was estimated using 
the iLandv & Szala'vl ( ll993l) estimator 

^ Nr{Nr - 1) GG(rp,n) Nr - 1 GRjrp, n) 
^ NciNc-l) RR{rp,7T) Nc RR{rp,n)^ 

where Nc is the mean galaxy density (or, equivalently, the to- 
tal number of objects) in the survey; Nr is the mean density of 
a catalog of random points distributed within the same survey 
volume; GGi/p, n) is the number of independent galaxy-galaxy 
pairs with separation between rp and rp + drp perpendicular to 
the line-of-sight and between n and n + dn along the line of 
sight; RR{rp,ji) is the number of independent random-random 
pairs within the same interval of separations and GR(rp, tt) rep- 
resents the number of galaxy-random cross pairs. 

The real-space correlation function ^(r) represents the ex- 
cess probabiht y for a given p air of galaxies to be observed at 
a separation r llPeebleslll980l) . It can be derived from ^(r^, n) 



using the formalism of lDavis'& PeehlesI lIlQSlh . computmg 
Wp(rp), the projection of ^(r^, n) along the line of sight. 

Wp(rp) = 2 ^(rp,n)dy^2 d(rl+y^)"^dy (2) 
Jo Jo 

In practice, the upper integration limit must be chosen finite, 
as to include the real signal, without adding extra noise which 
is dominant above a cer tain n. Based on d irect tests with the 
GallCS mock samples ( Blaizot et alJ200 5^. we choose a value 
^mnx = 20 h^'Mpc. This is high enough to sum the real sig- 
nal, without adding too much noise. It is compatible with the 
VVDS velocity measurement error of ~ 275 km.s"' . Assuming 
a power-law form for the real-space correlation function ^(r) = 
(I'/foY, the integral in equation |2] can be solved analytically 
in terms of Gamma functions, yielding a Wp{rp) which is it- 
self a pow er law. The fitting method of the measured Wp{rp), 
detailed in IPollo etalJil2nn.4 . provides the values of the corre- 
lation length ro and the slope y with their associated errors. A 
specific treatment have been applied to compute the clustering 
properties per galaxy types, as described in Sections 13.21 and 
l33l 

3.2. Biases and weights for each spectral class 

In order to correct for the observational biases introduced by 
the VVDS observing strategy in the VVD S-02h field, we use d 
the same correction scheme explained in IPoUo et al] ( l2005l) . 
The main biases are the incompleteness in redshift measure- 
ment and the complex geometry of the field. This results in 
a spatial sampling that changes as a function of position in 
the field. To account for galaxies with unknown redshift (ei- 
ther unobserved or discarded because of a low redshift quality 
flag), the observed galaxies are given a weight. The goal of this 
weight is to recover the true number of galaxy pairs at a given 
angular scale. 

The computation of this local weight takes into account the 
local number of galaxies (with and without a redshift) around 
a given galaxy in the spectroscopic sample. The parent photo- 
metric sample used for this computation must have the same 
properties as the spectroscopic one. For instance, the weight 
computed for the early type galaxies from the spectroscopic 
sample must make a reference only to the early type galax- 
ies in the underlying complete photometric sample, without in- 
cluding other galaxy types. In this approach, the parent pho- 
tometric catalogue used to compute the weight must therefore 
only contain galaxies with the same properties as those of the 
spectroscopic sample. 

Photometric redshifts of all the galaxies in the VVDS- 
02h have been derived from BV RI photometry dLe Fevre et alJ 
"20041 iMcCracken et Zl l2003l) . and ugriz band s from the 
CFHT L egacy Survey. Details are presented in lllbert et all 
din prep J) . The photometric redshifts and rest-frame magni- 
tudes were then used to give a spectral type to all these galaxies, 
according to the method explained in Section fT/H The com- 
plete parent photometric catalogue was then divided in four 
catalogues, one for each type. Note that the photometric red- 
shifts are only used to define a spectral type to each galaxy in 
the whole photometric catalog in order to compute the weights 
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associated to each galaxy type, in addition to the spectroscopi- 
cally measured galaxies. Photometric redshifts are not directly 
used here in the computation of the correlation function. 



3.3. Error measurement 

Even if our spectroscopic sample is the largest to date at 
this depth, and despite the angular size of the VVDS-02h 
(1750 arcmin^), the uncertainties associated to the measure- 
ment of the power-law parameters of the real-space correlation 
function ^(r) are largely dominated by the cosmic variance. 

To assess this, we have constructed 50 VVDS-02h mock 
surveys from the GallCS simulation jBlaizot et alJ l2005l) . 
GallCS is a semi-analytic model for galaxy evolution post- 
processed to a large cosmological N-body simulation. These 
mock surveys are quasi-independent and are built applying all 
the VVDS-02h observing biases: the complex geometry of the 
field, the spectroscopic selection function, the redshift mea- 
surement accuracy, the incompleteness (see lPoUo et al.l (l2005h 
for more information). We use the clustering variance in the 
mock samples to quantify the cosmic variance expected on the 
correlation length ro and the slope y of the real-space correla- 
tion function. 

Our method to assign a spectral type to a given galaxy does 
not work so well with the GallCS simulations. We were not 
able to recover a realistic redshift distribution for the four spec- 
tral types using the CFH12k BVRI apparent magnitudes pro- 
duced by GallCS. For example, there were no Type 1 (ellip- 
tical) galaxies beyond z~0.5 in most of the 50 GallCS cones. 
Conversely, the galaxy type classification given in for GallCS 
data do not match ours. This "Morphological" type classifica- 
tion is based on rest-frame B-band (bulge+burst)-to-disc lu- 
minosity ratio. This provides a huge "spiral" class (on aver- 
age ~90% of the full dataset), and two equivalent (in number 
of objects) "elliptical" and "lenticular" classes. These two last 
classes could match our type 1 class. But then, this is quite dif- 
ficult to divide the huge "spiral" class in our 3 later types. 

Computing error measurements from a bootstrap resam- 
pling method slightly under-estimates the real variance: if we 

consider the full VVDS-02h data set jLe Fevre et al 2005b) . 

the relative errors computed from GallCS are ~ 1 .5 - 2 times 
greater than error computed from bootstrap resampling. We ex- 
pect similar factor when dividing our sample per spectral type. 

Therefore the approach that has been chosen to compute 
measurement errors consists in taking randomly, in each of the 
50 mock VVDS-02h surveys, a proportion of galaxies equal to 
that found in the VVDS-02h spectroscopic sample for a given 
spectral class. Even if the clustering measured is the one of the 
"full" population with higher incertainties, the statistic of each 
"type" is representative of the real one. The variance from field 
to field is artificially increased by a broader range of cluster- 
ing properties in the simulated data, and therefore, this method 
is very conservative in the way that it is expected to produce 
errors for the observed data larger than in reality. 




0.1 1 10 1 a 3 4 5 6 7 8 



hrp (Mpc) hr^ (Mpc) 

Fig. 4. (left) Correlation function Wp(rp) as a function of red- 
shift for early-type galaxies with types 1 and 2. (right) Error 
contours (68%, 90% and 95% likelihood levels) associated to 
the measurement of the correlation length ro and slope y of the 
correlation function. 

4. Clustering evolution for each galaxy type 

4.1. Measurement of ro andy 

We have computed the correlation function ^(r^, n) and its pro- 
jection Wp(rp) in the VVDS-02h field, for each galaxy type of 
our classification, in increasing redshift slices. In order to im- 
prove the clustering signal, and depending on the redshift bin 
considered, we merged the two earliest-type galaxy samples Tl 
and T2 and the two latest types T3 and T4. Merging types 1 and 
2 would clearly dilute the clustering signal of type 1 galaxies 
if they are more strongly clustered than type 2. The measure- 
ments of Wp{rp) are presented in Figures|3to0 

The projection Wp{rp) has been fitted by a power-law on 
scales 0.1 < < 10 h 'Mpc. The values of the correlation 
length ro and the slope y of the real-space correlation func- 
tion are summarized in Table|2] A clear difference in clustering 
strength is observed between the 3 classes of galaxies: T1&2, 
T3 and T4. We find that on average the clustering strength 
of early type galaxies is 1.8 times larger than the clustering 
strength of type 4 galaxies. This is observed over the full red- 
shift range 0.2 < z < 1.2 and this result is significant. The 
probability that the correlation lengths of these two popula- 
tions are the same is indeed only 2.3% for the redshift bin 
z=[0.2,0.5], 2.7% for z=[0.5,0.7], 2.2% for z=[0.7,0.9], and 
14.5% for z=[0.9,1.2]. 

The evolution of the clustering length for the different 
classes of galaxies is presented in Figure|S] An overall increase 
in the clustering of all galaxy types is observed up to z ^ 1. 
In the range z=[0.9,1.2], there is a slight indication for a weak- 
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Fig. 5. Same as Figure|4]for galaxies with type 3. 



Fig. 7. Same as Figure|4lfor galaxies with type 3 and 4. 
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1.5 
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Fig. 6. Same as Figure|4]for galaxies with type 4. 



ening in the clustering of early type galaxies. The clustering 
length of late type / star-forming galaxies continues to rise up 
to redshift ~ 2. 

4.2. Relative bias and its evolution 

Since it is not possible to measure directly the real-space cor- 
relation function of the mass, any attempt to measure the bias 
b from galaxy samples has either to assume the knowledge of 



the background cosmological parameters (e.g. iMarinoni et alJ 
120051) . or to obtain the mass variance from other observables 
like the Cosmic Microwave Background (e.g. iLahav et"al] 
2002). 

The simplest definition of the bias b that links linearly the 
mass density field p of the mass to the galaxy density field n is 
(for a fixed scale R over which p and n are measured) 



6n 6p 



(3) 



where b is simply a number, independent on scale and on the 
local value of 6pl (p) hypothesis. Similarly, but in a statistical 
rather than deterministic way, the bias can be defined in terms 
of nns values of the galaxy and density fields 



^] ^bM 



(4) 



With this definition, b is simply a ratio if the square root of the 
variance of the galaxy and mass density fields measured over a 
given scale R, i.e. cr^(galaxies) - b (Tr (mass). In our case, we 
can measure directly the relative bias between two populations 
A and B, as 



bA 
be 



(5) 



where erg can be computed from the measured correlation func- 
tion as 



(6) 



which, for power-law form of ^(r), can be written as JPeeblesI 
II980) 



'"0 



8 Mpc/h 



(7) 
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- □ Types l&c2 
r# Types 3&4 




1 - 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

redshift z 

Fig. 8. Correlation length ro as a function of redshift for early (T1&2) and late (T3&4) type galaxies. Here the slope y is a free 
parameter. 



with 



' (3-y)(4-y)(6-y)2y ' 

Figure |9] shows the relative bias between the populations 
T1+T2 and T4 at different redshifts, measured using eq. 
in spheres with radius r = 8 h 'Mpc. The average value is 
b - 1.65 + 0.3 consistent with being constant up to the highest 
redshift investigated. This result seems to be mostly insensitive 
to the way the two populations are selected. For example we 
find a relative bias b = 1 .6 ± 0.3 when the red and blue samples 
are selected according to color bimodality (see next section). 
Note, however that by considering the intermediate-type popu- 
lation in our comparison, we find a relative bias b = 1.3 + 0.2 
between T1-I-T2 and T3, which is marginally (~ Icr) lower than 
between T1-I-T2 and T4. 

This is in agreement with what is found locally for nearly 
the sa me color-selected populations. For example Wil d et alJ 
l l2005h using the 2dFGRS find on the same 8 h 'Mpc scale, 
b ~ 1.5. Moreover, at higher redshifts, this result is in agree- 
ment with what has been found in a previous VVDS study 



jMarinoni et aljEoO^ . from the same data but using a differ- 
ent method. 



5. Rest-frame color bimodality 

It has become clear in recent years t hat the galaxy rest - 
frame colors show a bimodal distribution dStrateva et al.l200ll) . 
i.e. that in a (rest-frame) color-magnitude diagram galaxies 
tend to segregate between a "red sequence" (similar, but less 
tight than that observed for cluster galaxies) and a "blue 
cloud". This behavior seems also to be present at high redshift 
(Bell et al." 2004:"Gi allongo et alJ l200.5Y From the SDSS data, 
.Strateva et al... C2001il) find that the blue cloud contains mainly 
late (spiral) morphological types while the bulk of the red se- 
que nce consists of bulge-dominated, early-type galaxies (see 
also lWeiner et aL.2005.) . This provides a natural recipe to split 
galaxies into two populations of "red" and "blue" objects and 
there have been already attempts to study the clustering of 
galaxies separately for the two cl asses define d in this way at 
intermediate redshifts (e.g. Phleps et alJ2005|) . 
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Table 2. Measurement of the clustering length for galaxy subsamples ranging from z-0.2 to z-2, letting the slope y free or fixing 
the slope to the average value measured for a given type. The absolute Mg ^^ magnitude and number of objects are indicated. 



type 


number of 


redsliift 


effective 




5log(h) 


'■o(z) 


r 


ro(z) at fixed y 




galaxies 


range 


redshift 


mean 


median 


(h-'Mpc) 




(h-'Mpc) 


1 


164 


[0.2-0.6] 


0.427 


-18.799 


-19.020 


3 14+0.76 

-^•^ -0.94 


7 41+0.33 
^•^^-0.37 


3.51!U;^^ (7 = 2.15) 




332 

237 


[0.6-1.0] 
[0.8-1.2] 


0.779 
0.964 


-19.755 
-20.324 


-19.819 
-20.347 


4 ^S+0.68 

3 89+' "'* 


2 01+° "* 

Q29 

9 04+0.80 


T QT+0.55 

-^•'-^-0.74 

3.62+°f 


2 


299 


[0.2-0.6] 


0.409 


-18.352 


-18.419 


9 7Q+0.42 
^- '="-0.44 


1.88-- 


2.76!;j§(y= 1.90) 




420 
358 


[0.6-1.0] 
[0.8-1.2] 


0.796 
0.972 


-19.525 
-20.059 


-19.503 
-20.115 


4 Q9+0.72 
4 Q 5+0.96 


1 OO+0.27 
i-^^-0.20 

2 00+°^** 

^•""-0.26 


4 00+0.58 
^■°^-0.58 

4 62+° '*^ 



1&2 



437 
470 
473 
439 
191 



507 
598 
534 
591 

289 



307 
373 
390 

371 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 
[1.2-2.0] 



0.346 
0.606 
0.810 
1.027 
1.394 



-17.478 
-18.665 
-19.296 
-19.843 
-20.756 



-17.332 
-18.601 
-19.212 
-19.835 
-20.751 



2 69+"-^" 

^•"'-0.54 
Q 71+0.60 
-'•^^-0.69 
Q 74+0.62 
-'•^^-0.65 
O QC+0.44 
^.SS_o44 

3.60+?-^^ 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 
[1.2-2.0] 



0.366 
0.604 
0.807 
1.036 
1.394 



-17.255 
-18.277 
-18.944 
-19.553 
-20.441 



-17.170 
-18.136 
-18.829 
-19.494 
-20.397 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 



0.346 
0.604 
0.814 
1.032 



-18.327 
-18.996 
-19.719 
-20.333 



-18.512 
-19.020 
-19.635 
-20.383 



1 7 1 +0.30 
^- ' ^-0.31 
7 Q 1 +0.40 
^•-'^-0.40 

2 46+°-'" 

9 Qc+0.36 

Q Qn+0.64 
-^•--^ -0.87 
Q 9 1 +0.58 
■-^•^^-0.60 
4 1 1 +0.66 
^•^^-0.76 

4.26+!;-,^I 



1.83-- 

I Q 1+0.24 
^•'^-0.15 

1 56+°-22 

1 QQ+0.36 

1 7Q+0.35 
^•'•'-046 



2-73!U;:^(y= 1.78) 



n TT+0.54 

J.JJ Q^g 



3.21+ 



T80^ 
1.76!°;!^ 

I Q9+0.31 
^•'^-0.19 
1 Q4+0.I8 

1 Q9+0.36 
^•^ -0.32 
1 Q4+0.27 
-0.29 
1 Q4+0.I9 
^•'^-0.12 

1 76+° ""^ 

I Q9+0.31 



I +0.60 
-0.60 
9 89+0.49 

3 65+' °2 

-^•"-^-3.04 



9 9Q+0.42 
2.28 1,3, 

7 40+0 53 
^•^^-0.58 
7 Q 1+0.40 
^•^^-0.37 
T 97+1.01 
-'•^'-2.35 



1.86) 



3-21!!!::!! (7 = 1. 

4 14+0.53 
^•^ -0.57 
4 41+0.62 
^•^^-0.69 
4 ^9+0.79 



3&4 


943 


[0.2-0.5] 


0.357 


-17.358 


-17.230 


9 1 9+0.35 
^•^^-0.36 


1 70+0. 14 
i-'*-0.09 


2.12!«|(y= 1.76) 




1068 


[0.5-0.7] 


0.605 


-18.448 


-18.321 


2 66+°-''^ 

^•■^"-0.45 


1 74+0- '5 

^ - ' -0.09 


2 66+° '" 

^•""-041 




1007 


[0.7-0.9] 


0.809 


-19.109 


-19.008 


9 40+0.42 
^•^'-0.42 


1 67+°'^ 
^■"'-0.11 


2.46!°;^^ 




1030 


[0.9-1.2] 


1.032 


-19.677 


-19.626 


9 5O+0.25 
z..jo_Q22 


1 sA+o.ii 

^•^°-0.08 


2.58!°-5,« 




480 


[1.2-2.0] 


1.394 


-20.566 


-20.510 


Q 49+0.66 
■-^ -0.73 


1 QQ+0.24 
1.77 


T 94+0.90 

-^•^ -1.23 



Table 3. Relative bias between early-type and late-type galaxies (from spectral type classification), and "red" vs. "blue" galaxies 
(separated using the color bi-modahty) 



redshift relative 



range 



bias 

1 1O+0.23 

1 97+0.26 
^■■^'-0.28 
I OQ+0.26 
^■-^'-0.27 

1 45+°^" 



T1&2/T3 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 



T1&2/T4 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 



1 76+°" 

1 . / U_Q ^-J 

1 69+°-^° 

1 0^+0.43 

^•^-'-0.47 
1 47+0.37 



1 69+"« 
1 63+°'^ 

l.UJ 1,34 

1 45+°2« 

l.tj Q3g 

1 45+°" 



"red"/"blue" 



[0.2-0.5] 
[0.5-0.7] 
[0.7-0.9] 
[0.9-1.2] 



It is therefore interesting to check how our analysis as 
a function of SED-defined types compares to a simpler sub- 
division following the bimodal distribution. 

We have thus split our data into a blue and red sample, 
fo llowing the rest-frame color-magnitude relations suggested 
bv lGiallongo et alJ ( l2005l) and measured the projected function 
Wp(rp) and the best-fit power-law correlation function ^(r) in 
different redshift bins. The resulting evolution of the correla- 
tion length is shown in Fig.^| with the measured values for 
and y summarized in Table [S] The associated errors have been 
computed in the usual way, constructing blue and red mock 



catalogues from the GallCS simulations, keeping the same pro- 
portion of red (blue) objects as in the real VVDS data. 

Similarly to what we found for early-type classes, red- 
sequence galaxies exhibit a larger clustering length (ro ~ 
4 h 'Mpc), with little dependence on redshift. Comparison to 
Fig. |S] indicates that the classification into rest-frame red and 
blue galaxies is substantially equivalent to our classification 
into spectral late-types, producing samples with similar clus- 
tering evolution properties. 

The relative bias between red and blue galaxies is measured 
to be on average b = 1.6 + 0.3 between z ~ 0.3 and z ~ I 
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Fig. 9. Relative bias between populations of galaxies for a Fig. 10. Correlation length ro as a function of redshift for the 
8 h^'Mpc scale: ^a/^b = ^l^t- h\vLt samples. Here the slope y is a free parameter 



(see Tabled, a value comparable with that derived using the 
spectral type classification. 

In summary, using the bi-modal split of the population 
we find results consistent with our analysis per spectral types. 
However, we note that the simple rest-frame color-magnitude 
cut is a less precise way of selecting truly early-type galaxies, 
at least for z > 0.7, since the full multi-color information is 
used to provide the spectral type classification. 

6. Discussion and conclusions 

The fundamental result of this investigation is to have clearly 
established that at least up to z ^ 1.2 early-type galax- 
ies continue to be more strongly clustered than late-type 
galaxies. This extends to h igher redshifts earher results ob- 
tained in the local Universe (iLovedav Tresse & Maddoxll999t 
iNorberg et al.'2002': 'Z ehavi et al. "2002'. 'lOOSY at redshifts z ^ 

0. 6 tPhleps et al.-2005<). and auantitativelv confirms early indi- 
cations atz ^ 1 (iPhleps & M eisenheimer 2003). 

Our specific results can be summarized as follows: 

1 . The clustering of sub-L* early-type galaxies is stronger than 
the clustering of late-type galaxies up to z ^ 1 .2; a similar be- 
havior is obtained for "red-sequence" and "blue-cloud" galax- 
ies selected using rest-frame colors. 

2. Luminous, early-type galaxies are already strongly clustered 
at z ~ 1 with 3.5< ro < 5 h^'Mpc. 

3. The clustering of late-type galaxies is low at all epochs with 
ro ^ 2.1 -3.5 h-'Mpc. 

4. The overall clustering of bright late-type galaxies at z - 1 .4 
is relatively strong with ro = 3.4 h^'Mpc 

5. The bias between the early-type and late type population 



is nearly constant, beariylbiate - 1-7 over the redshift range 
0.2<z<1.2. 

A physical interpretation of these results and comparison to 
local values require to keep in mind that we are working with a 
flux-limited sample, and thus that we are measuring galaxies of 
diff'erent intrinsic luminosity at different redshifts. At z ~ 1 the 
VVDS targets early-type galaxies with -22 < Mb{AB) < -19 
and an average Mb(AB) = -20.3 (see Figure |3}- Early-type 
galaxies observed in the 2dFGRS with similar luminosities 
have a n observed rp ^ 6 - 7.5 h^'Mpc, increasing with lumi- 
nosity dNorberg et alJl2002h . The clustering amplitude of ^(r) 
for early-type galaxies has therefore evolved by no more than 
70% between redshift ~ 1 and present. 

To understand the implications of this observation, we have 
compared our observed evolution of ro to the predicted evo- 
lution of the clustering length for the dark matter in a A cos- 
mology ( using the expe cted evolution of the variance as de- 
scribed in Peebles (1980)). As shown in Figure^J if one as- 
sumes that the early-type population is simply tracing in an 
unbiased way the growth of clustering in the mass (as would 
have been the case if no evolution in number nor in lumi- 
nosity occurred) their clustering would increase by a factor of 
cTaiz = l)/cr8(z = 0) ^ 0.6, which is close to what we see 
in Figure This result is in agreement with the observation 
of a marginal clustering evolution of early -type galaxies since 
z ~ 1.5 (Cima tti et al.ll2004t ISaracco et alj E004). i.e. with flie 
idea that the bulk of the spheroid population was already in 
place at these early redshifts. If a significant fraction of these 
galaxies was assembled and added below z ~ 1, we would ob- 
serve not only an increase in their number density, but also a 
weakening of their clustering signal with respect to the under- 
lying mass: it is natural to think that, forming typically in lower 
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Table 4. Number of galaxies and mean B rest-frame absolute magnitude (h= 1 ) of the sample separated using the color bimodality. 
The number of "red" and "blue" galaxies classified as early (T1&2) and late (T3&4) spectral types are identified. We note that 
the populations of "red" and T1&2 as well as the populations of "blue" and T3&4 galaxies are mostly the same, except for a 
significant cross-overpopulation which could explain the observed difference in clustering (see text). 



redshift 


red 




blue 




range 


T1&2 


T3&4 


T1&2 


T3&4 


0.2 - 0.5 


221 


49 


86 


895 




-18.181 


16.403 


-18.701 - 


17.411 


0.5 - 0.7 


253 


24 


121 


1044 




-18.892 


17.646 


-19.215 - 


18.467 


0.7 - 0.9 


313 


42 


77 


965 




-19.713 


18.468 


-19.748 - 


19.138 


0.9- 1.2 


283 


13 


88 


1016 




-20.263 


19.314 


-20.557 - 


19.679 


1.2-2.0 


41 





94 


480 




-20.926 





-21.078 - 


20.567 



Table 5. Measurements of the clustering length and slope for galaxy subsamples divided following the color bimodality 



sample number of redshift effective Ms,g - 5!og(h) ro(z) y 





galaxies 


range 


redshift 


mean 


(h-'Mpc) 




red 


270 
277 


0.2 - 0.5 
0.5 - 0.7 


0.345 
0.610 


-17.858 
-18.784 


3 84+"-»'' 

A TC+0.65 


1.88!-^ 

1 Qq+0.22 

^•''^-o.is 




355 
296 


0.7 - 0.9 
0.9- 1.2 


0.809 
1.022 


-19.565 
-20.221 


T 78+0.70 
^•'°-0.74 
Q 79+0.63 
' -0.67 


1 87+0.28 
1.8/_D22 
1 eQ+0.32 


blue 


980 
1164 
1042 
1105 
574 


0.2 - 0.5 
0.5 - 0.7 
0.7 - 0.9 
0.9- 1.2 
1.2-2.0 


0.357 
0.604 
0.811 
1.035 
1.404 


-17.521 
-18.544 
-19.183 
-19.749 
-20.651 


z.O'_0 35 
7 cc+0.39 

7 1 O+0.32 
2.18 Q 

7 4Q+0.28 
^•^'-0.22 
Q 74+0.63 
-0.70 


1 67+"''' 

i-°°-0.09 

1 40+""' 

1 84+"'" 
^■°^-o.io 

1 Q 1+0.20 
^•'^-0.18 



density peaks than their high-redshift companions they would 
in fact be less massive and less clustered. 

The lower clustering level observed for early-type galax- 
ies at z < 0.7, is best interpreted as a luminosity effect, given 
the faint typical luminosity of galaxies in the first bins (e.g. 
{Mb(AB)} = -18.3 atz = [0.2,0.5]). 

The trend observed for blue star-forming galaxies of simi- 
lar luminosities is also very interesting. Below redshift z - 1, 
their clustering is observed to be weak and continuously de- 
creasing from high redshift to the present. As shown in the 
two right-top panels of Figure ^2 however, this is consis- 
tent with simple gravitational growth: the clustering amplitude 
measured at z ~ 1 for these galaxies is coherent with the evo- 
lution of the same structure traced by the same galaxies at 
z ~ 0. The point at z = 1.4 may indicate that star forma- 
tion activity at z > 1 is found in progressively more strongly 
clustered galaxies than at present, i.e. plausibly in more mas- 
sive galaxies in higher density peaks. In our sample, galaxies 
at z = [1-2,2] have a mean Mb{AB) — -20.5 and a cluster- 
ing length ro = 3.5 h 'Mpc. Taking into accou nt 1.5 magni- 
tudes of B-band brightening ( IZuccaetalJl2005l) these galax- 
ies would have today Mb(AB) — -19. Late-type galaxies with 
similar luminosities in the Stromlo-APM, 2dFGRS or SPSS 
have rn - 2.9 - 4.2 h~'Mpc ('Lovedav Tresse & Maddox 1999; 
iNorberg etai]l2002t l Zehavi et al. .2002.) . The clustering am- 
plitude of late-type star-forming galaxies with comparable lu- 



minosities, therefore remains roughly constant since z ~ 1.5 
(Figure [m. A straightforward explanation to produce the ob- 
served result is that as redshifts get smaller, star formation is 
shifting to smaller-mass, less clustered galaxies. 

It is difficult to compare correlation properties of samples 
with different mean/median magnitudes because of the well 
know n dependence of clustering on luminosity (Iguzzo et"al 
'1997| iBudavari et all l2003t IZehavi et alJ l2002t IPollo etal 
2006;). Nevertheless, we can note that comparing our values 
with other clustering measurements for type- or color- selected 
galaxies at redshifts beyond 1.5, we find for example that star- 
forming LBGs at z ~ 3 ( Adelberger et al. 2004) are even more 
clustered (ro ~ 4.5 h^'Mpc) than our star-forming galaxies at 
z ~ 1.5. This is consistent with the VVDS results at lower red- 
shifts. One can speculate that star formation is moving to higher 
and higher density environments, when moving back in time. 
Similarly, red selected samples of EROs at z ~ 1 .5 or red K- 
selected galaxies with 2 < Zphot < 4 are significantly more clus- 
tered than the e arly type galaxies at z - 1 in the VVDS with ro 
up to 8 h 'Mpc jDaddi et alJ2002ll2003i) . These galaxies seem 
to be populating even more extreme density peaks than those 
we have been studying here at z ~ 1 . 
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Fig. 11. Clustering length evolution for early (left panel) and late (right panel) spectral type galaxies. Measurements from the 
VVDS at high redshifts (open symbols) are compared to the local measurements obtained for galaxies in the same absolute 
luminosity ranges in the 2dFGRS (filled symbols) surve ys (J^orberg et al..200Z) . with an evolution correction of 0.5 magnitudes 
for early types and 1 magnitude for late types at z ~ 1 jZucca et alJl2005h . The dashed curves show the predicted gravitational 
growth of the dark matter, with amplitude normalized to erg = 0.9 at the current epoch. 
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